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We derive the equations of motion in metric-affine gravity by making use of the conservation laws 
obtained from Noether's theorem. The results are given in the form of propagation equations for the 
multipole decomposition of the matter sources in metric-affine gravity, i.e., the canonical energy- 
momentum current and the hypermomentum current. In particular, the propagation equations 
allow for a derivation of the equations of motion of test particles in this generalized gravity theory, 
and allow for direct identification of the couplings between the matter currents and the gauge 
gravitational field strengths of the theory, namely, the curvature, the torsion, and the nonmetricity. 
We demonstrate that the possible non-Riemannian spacetime geometry can only be detected with 
the help of the test bodies that are formed of matter with microstructure. Ordinary gravitating 
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spacetime geometry. Thereby, we generalize previous results of general relativity and Poincare gauge 
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I. INTRODUCTION 



The relation between the field equations and the equations of motion within nonlinear gravitational theories has 
been subject to many works. The intimate link between these equations is one of the features of general relativity 
which distinguishes it from many other physical theories. The fact that, in contrast to linear field theories, the 
equations of motion need not to be postulated separately, but can be derived from the field equations, has been 
investigated shortly after the proposal of the theory. From a conceptual standpoint the derivability of the equations 
of motion is a very satisfactory result, since it reduces the number of additional assumptions in the theory. 1 The 
earliest accounts of this feature of general relativity can be found in the works of Weyl 0, Eddington [H, as well 
as Einstein and Grommer Nowadays this is customarily addressed as the "problem of motion" in the context of 
general relativity and other nonlinear field theories. 2 

One may distinguish between two conceptually different methods. Both were employed in the derivation of the 
equations of motion within the theory of general relativity. One of them goes back to the works of Einstein et al. 
[H, [|| and is based on the vacuum field equations of the theory. Within this method matter is modeled in the form of 
singularities of the field and only the exterior of bodies is considered. The second method, usually attributed to Fock 



'Electronic address: dirk.puetzfeld@astro.uio.no URL: http://www.thp.uni-koeln.de/~dp 
^Electronic address: yo@thp.uni-koeln.de 

* Also at Department of Theoretical Physics, Moscow State University, 117234 Moscow, Russia 
1 The following german quotes are taken from []J (translation by the authors): 

• "[. . .] Es sieht daher so aus, wie wenn die allgemeine Relativitatstheorie jenen argerlichen Dualismus bereits siegreich iiberwunden 
hatte. [. . .]" , 

"[. . .] It looks like the general theory of relativity has victoriously overcome this annoying dualism. [. . .]". 

• "[. . .] Der hier erzielte Fortschritt liegt aber darin, dafi zum ersten Male gezeigt ist, daft eine Feldtheorie eine Theorie des 
mechanischen Verhaltens von Diskontinuitdten in sich enthalten kann. [. . .] " , 

"[. . .] The progress achieved in this work is that for the first time we have shown that a field theory can contain the theory of the 
mechanical behavior of discontinuities. [. . .]". 



2 A historical account of works can also be found in 0, U 0, . 
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1 1 01 ] , makes use of the differential conservation laws of the theory and also allows for a consideration of the interior of 
material bodies. In this work we are going to utilize the latter method; i.e., we base our considerations on differential 
identities derived from the symmetry of the action via Noether's theorem. 

In addition, we make use of a multipolc decomposition of the matter currents. This allows for a systematic study 
of the coupling between the matter currents and field strengths of the theory at different orders of approximation. 
Multipolc methods have been intensively studied in the context of the problem of motion since the early work of 
Mathisson [Til ]. I n table HI we provide a corresponding chronological overview. 3 

TABLE I: Timeline of works which deal with the problem of motion and 
multipole approximation schemes. 



Year Reference 



Comment 



1923 Weyl [2] 

1927 Einstein and Grommer [1] 

Lanczos [131 
1931 Mathisson [14, JJ5, 16] 

1937 Robertson [17] 
Mathisson [11] 

1938 Einstein et al. [8, 9] 

1939 Fock [10] 

1940 Papapetrou [18] 

1941 Lanczos [19] 

1949 Infeld and Schild [20] 
1951 Papapetrou [21] 

Papapetrou |23] 
1953 Papapetrou [24] 

Goldberg [25J 
1955 Meister and Papapetrou [261 
1957 Infeld [27] 

1959 Kerr [28, 29] 

Fock [30J 
Tulczyjew [31] 

1960 Infeld and Plebanski [32] 
Kerr [33] 

Synge [34] 

1962 Goldberg [5] 

Havas and Goldberg [35] 
Tulczyjew and Tulczyjew [36] 
1964 Taub [37] 

Dixon [38] 
Havas [39] 

1969 Madore [40J 

1970 Dixon [41, 42] 
1973 Liebscher [43, 44] 



Mentions the link between the equations of motion (EOM) and the 
field equations. 

Show that the field equations contain the EOM in GR (for a special 
case) . 

Early investigation regarding the problem of motion, treated as 
boundary value problem. 

Systematic account of the problem of motion in GR, one of the first 
authors who makes use of the <5-function in this context. 
Test particle EOM from divergence condition. 

Possibly the earliest work utilizing a multipole method in the deriva- 
tion of the EOM. 

Derivation of the EOM outside of material bodies. 
Systematic slow motion approximation. 

Gravitational interaction of particles using the multipole method. 
Test particle EOM via Gaussian integral transformation. 
Derive the geodesic motion of test particles for empty space. 
EOM for pole-dipole test particles in GR (see also the later work 

HI). 

Derivation of the EOM utilizing a method in the spirit of |10( | . 
Review of the problem of motion in GR. 
Relationship of EOM and covariance of a field theory. 
EOM and coordinate conditions in GR. 

Review of approximation methods, derives EOM using Einstein- 
Infeld-Hoffmann (EIH) method, relaxes harmonic coordinate con- 
dition, <5-function as source. 

Systematic post-Minkowskian treatment I + II (fast motion 
approximation) . 

Systematic slow motion/weak field approximation. 

Test particle EOM via a simplified version of Mathisson's method. 

Review of the EIH method. 

Approximation of the quasistatic case, review of three approxima- 
tions schemes. 

Integral conservation laws, EOM for mass center, energy- 
momentum pseudotensor definition. 

Review of the problems connected with the EOM in GR and the 
EIH method. 

Derive single-pole EOM by using Mathisson's method. 
Covariant formulation of a multipole method in GR. 
Test particle EOM in a coordinate independent manner using Pa- 
papetrou's method. 

Covariant multipole method for extended test particles in GR. 
Generalized version of Mathisson's method in afline spaces. 
EOM for extended bodies using a multipole method which differs 
from the one of [2l| . 

Extended bodies within a multipole formalism. 

EOM for pole particles in non-Riemannian spaces using the method 

in EU, see also '451. 



3 An extended version of this table, also including works in the post-Newtonian and post-Minkowskian context, can be found in [T3l 
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1974 Papapetrou [46J 



Review of the derivation of the EOM of a single-pole test particle 
in GR. 

Review of the multipole formalism in GR in the context of extended 
bodies. 

Generalization of the Papapatrou equations to Poincare gauge 
theory. 

Multipole method for the derivation of the EOM for extended 

bodies. 

EOM review. 



1979 Dixon [47] 



1980 Yasskin and Stoeger [48] 



Bailey and Israel [49] 



1987 Damour [7] 



In this paper, we work out the equations of motion within a multipole formalism for a generalized gravitational 
theory known as metric-afHne gravity (MAG) [50]. In the theory of general relativity, the mass, or more precisely the 
energy-momentum, of matter is the only physical source of the gravitational field. The energy-momentum current 
corresponds (via the Noether theorem) to the local translational, or the diffeomorphism, spacetime symmetry. In 
MAG, this symmetry is extended to the local affine group that is a semidirect product of translations times the local 
linear spacetime symmetry group. Correspondingly, there are additional conserved currents describing microscopic 
characteristics of matter that arise as physical sources of the gravitational field. In continuum mechanics [HTL l52l . 
[531 [5^ , I551 . HH , such matter is described as a medium with microstructure. In physical terms this means that the 
elements of a material continuum have internal degrees of freedom such as spin, dilation, and shear. The three latter 
microscopic sources are represented in MAG by the irreducible parts (that correspond to the Lorentz, dilational 
and shear-deformational subgroups of the general linear group) of the hypermomentum current. Fluid models with 
microstructure were extensively studied within different gravity theories (including MAG), see, e.g., [57l l58l l59l l6fl . 

The metric-affine theory naturally generalizes the Poincare gravity theory |62l . |63| in which the mass (energy- 
momentum) and spin are the sources of the gravitational field. The geometry that arises on the spacetime manifold is 
non-Ricmannian, it is known as the Riemann-Cartan geometry with curvature and torsion. In MAG, this geometrical 
structure is further extended to the metric-affine spacetime with curvature, torsion, and nonmetricity. The resulting 
general scheme of MAG embeds not only Poincare gravity, but also a wide spectrum of gauge gravitational models 
based on the conformal, Weyl, de Sitter, and other spacetime symmetry groups (for an overview, see 50], for example). 
This fact makes the analysis of the equations of motion in MAG especially interesting, with possible direct physical 
applications for all the gravitational models mentioned. 

The energy-momentum current and the hypermomentum current (spin + dilaton + shear charge) are the sources 
of the gravitational field in MAG. Accordingly, test bodies that are formed of matter with microstructure have 
two kinds of physical properties which determine their dynamics in a curved spacetime. The properties of the first 
type have microscopic origin; they arise directly from the fact that the elements of a medium have internal degrees of 
freedom (microstructure) . The properties of the second type are essentially macroscopic; they arise from the collective 
dynamics of matter elements characterized by mass (energy) and momentum. More exact definitions will be given 
later, but the qualitative picture is as follows. The averaging of the microscopic hypermomentum current yields the 
integrated spin, dilaton, and shear charge of a test body. In addition, the averaging of the energy-momentum and of 
its multipole moments gives rise to the orbital integrated momenta. In Poincare gravity, there is only one relevant 
first moment, namely, the orbital angular momentum. It describes the behavior of a test particle as a rigid body, 
i.e., its rotation. In metric-affine gravity, one finds, in addition, the orbital moments that describe deformations of 
body. These are the orbital dilation momentum (that describes isotropic volume expansion) and the orbital shear 
momentum (that determines the anisotropic deformations with fixed volume). The three together (orbital angular 
momentum, orbital dilation momentum, and orbital shear momentum) comprise the generalized integrated orbital 
momentum. In this paper, we compare the gravitational interaction of the integrated hypermomentum to that of the 
integrated orbital momentum of a rotating and deformable test body. Thereby, we generalize the previous analysis 
[48| in which the effects of the integrated spin were compared to the effects of the orbital angular momentum of a 
rotating rigid test body. 

The paper is organized as follows. In section [TT] we recall some basic facts about the gravity theory under consid- 
eration, namely, metric-affine gravity. This is followed by a discussion of the conservation laws within this theory in 
section Mil which form the basis for the derivation of the equations of motion. We then work out the explicit form 
of the propagation equations in sections IIVI and [V] In section IVII we provide some relations between the different 
definitions of momenta within the multipole formalism. We discuss our findings in section fVIII and present an outlook 
on the open questions within this field. Our notation and conventions are summarized in appendix [A] A table with 
the dimensions of all quantities appearing throughout the work can be found in appendix [Bl 
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II. METRIC- AFFINE GRAVITY 



Metric-affine gravity represents a gauge-theoretical formulation of a theory of gravitation which is based on the 
general affine group A(A, R), i.e., the semidirect product of the four-dimensional translation group i? 4 and the general 
linear group GL(4,R). For a review of the theory see HU, and references therein. In such a theory, besides 
the usual "weak" Newton-Einstein-type gravity, described by the metric of spacetime, additional "strong" gravity 
pieces will arise that are supposed to be mediated by additional degrees of freedom related to the independent linear 
connection T a ^ . Alternatively, the strong gravity pieces can also be expressed in terms 4 of the nonmetricity Q a p and 
the torsion T a . The propagating modes related to the new degrees of freedom are expected to manifest themselves 
in the non-Riemannian pieces of the curvature R a @ . The existence of such modes certainly depends on the choice 
of the dynamical scheme, or in technical terms, on the choice of the Lagrangian. The simplest generalization of the 
linear Hilbert-Einstein Lagrangian leads to a model with contact interaction. However, quadratic Yang-Mills-type 
Lagrangians describe a wide spectrum of non-Riemannian propagating gravitational modes. This is revealed, for 
example, by studies of generalized gra vitational waves in models with torsion (65l . l66l . l67l . l68l l69l . l70l . l7l| and in models 
with torsion and nonmetricity [z3, M 0, [zl M, IzS H [zS HI . 

In a Lagrangian framework one usually considers the geometrical "potentials" (metric g a p, coframe 1-form d a , 
connection 1-form Tj 3 ) to be minimally coupled to matter fields, collectively called tp, such that the total Lagrangian, 
i.e., the geometrical and the matter part, is given by 

Ltot = L (g af3 , iT, Q a p, T a , Rj) + L mat (g af3 , 0°, V> Dip) . (1) 

Here D = d + l a p F Q ' 3 , with t a p denoting the generators of the linear transformations (namely, Si/j = e 13 a l a p ip, where 
e^a are the infinitesimal parameters). With the following general definitions for the gauge field momenta 

the field equations of metric-affine gravity take the form 

{S/Sg a p) DM af3 - = (3) 

(S/5$ a ) DH a - E a = E a> (4) 

(S/STJ) DH a p — E a f3 — A"^, (5) 

SL , . 

(matter) — = 0. (6) 

Sip 

On the right-hand side (rhs) of the field equations we have the physical sources: the metrical energy-momentum a a P , 
the canonical energy- momentum S a , and the canonical hypermomentum A a p currents of the matter fields 

a - =2 ^7' E «-=^' a ^-=^f- (7) 

On the left-hand side (lhs) there are typical Yang-Mills-like terms governing the gauge gravitational fields, and the 
corresponding terms that describe the currents of the gauge fields themselves that arise due to the nonlinearity of 
the theory. The metrical energy-momentum, the canonical energy-momentum, and the canonical hypermomentum 
currents of the gauge gravitational fields are introduced by 

mQ/3 ._ 9 9L _ 8L _ PL 

^? Ea -^ E ' 3 -wj- (8) 

MAG has a wide gauge symmetry group. With the help of the Noether theorems for the diffeomorphism symmetry and 
for the local linear symmetry, one can verify that [provided the matter field equations ((6J arc fulfilled] the following 
identities hold: 

E Q = eajXmat - (e a \Dip) A 9 q£T - (e a \ip) A ^" at , (9) 



4 Please see appendix fS] on page !22l for the definitions of the objects in this section and a short summary of our conventions. 
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E a = e a \ L + (e Q J T^) A Hp + (e Q J Rp 1 ) A 77 /3 7 + ~ (e Q J Af'' 7 , (10) 

= -$ a /\H p -M a p, (11) 

A^ = (f^)A^, (12) 

££ Q = (ejT' 3 ) AS /3 -i(e Q JQ^)a^ + (e ct Ji 1 , /3 7 )AA' 3 7 , (13) 

DA a ,g = g Pl a a ~< i3 a A E^. (14) 



The gauge symmetry and the corresponding Noether identities play an essential role in MAG. The most important 
result is as follows: It can be shown that, by means of (flTj)) - (114l) . the field equation ([3]) is redundant. It is a consequence 
of the two other MAG field equations (|4|) and ([5]) and of the Noether identities. The explanation is straightforward: 
One can use the local linear transformations of the frames to "gauge away" the metric g a p by making it equal to the 
constant Minkowski metric diag(l, —1, —1, —1) everywhere on the spacetime manifold. After doing this, equation ([3]) 
is trivially solved, and one needs to solve only the remaining equations ^ and ([5]) to determine the coframe "d a and 
connection T a @ . 

There are many nontrivial exact solutions for different MAG models ranging from black holes, gravitational waves, 
to cosmological models known in the literature. Nearly all of the corresponding references can be found in the works 

[II IE [sill. 



III. CONSERVATION LAWS 



An up-to-date discussion of the conservation laws within metric-affine gravity can be found in the recent work 
[§ij . In the following sections IIII AflTTCl we recall the conservation laws for the canonical energy- momentum and 
hypermomentum. These conservation laws serve as starting point for our subsequent derivation of the propagation 
equations for the multipole moments of the matter currents. In IIII CI we make contact with Poincare gauge theory, 
which represents the special case of metric-affine gravity for which the distorsion, i.e., the difference between the full 
and the metric-compatible connection, reduces to the antisymmetric contortion, and the hypermomentum reduces to 
the spin current. 



A. Energy-momentum conservation 



The Noether theorem for the diffeomorphism invariance of the matter action yields the conservation law of the 
energy-momentum current 

D (S Q - Af p ea\N^) = (e a \R/ - L a N^ A AV (15) 

{> ,{} (},,,, 
Here L ^ = £J D + DS,\ is the (Riemannian) covariant Lie derivative. 

After we substitute the components from (I A5|l - (IA10|1 . we finally find the tensor form of the conservation law (|T5)) : 



Vj (TV - N lkl A^') = ( { R llkl - vW iW )A fc y (16) 



This can be identically rewritten as 



where we denoted 



Vj Tj = Ryu A fcy + N lkl VjA^' , (17) 



{} {} {} 

Rijkl '■= Rijkl — ViNjkl + VjNikl- (18) 
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B. Hypermomentum conservation 



The Noether theorem for the local GL(4, i?)-invariance of MAG yields (on the mass shell, i.e., when the matter 
satisfies the field equations): 

DA a (} + d a A E - <T a a = 0. (19) 

Here the last term describes the metrical energy-momentum 4- form defined in equation ([7]). By the introduction of 
local coordinates for the corresponding components, 

a °- (i = t af3 n, (20) 

we can rewrite the Noether identity (|19|) in tensorial form: 

Vj A fey - N l3 k A jH + N jU A k i:j + T lk - t kl = 0. (21) 
Taking the antisymmetric part, we find: 

A [kl]: > = Nij [k A W]i + N j[k ^A l] ij + T [kl] . (22) 



C. Recovering Poincare gauge theory 



The case of the Poincare gauge theory is recovered when the difference of the connections reduces to the anti- 
symmetric contortion N a p = K a p = K\ a m, whereas the hypermomentum reduces to the antisymmetric spin current 
A a p = T a p — T[ a py 

With the help of (|2"2"]) , we then immediately find 



{} 



Kiki V jT hlj — Kiki T kl + (K in i Kj k n — Kj n i K lk n ) t 
Substituting this into (|17p . and rearranging the rhs, we have 

V, I, J 11,11 T kl i+K ikl T kl . 
Here the total Riemann-Cartan curvature is recovered in the first term on the rhs: 



klj 



{} {} {} 
Rijki — Rijki — ViKjki + V jKiki 



Kinl K 



jk 



^ jnl ik 



in complete agreement with (IA2[) . 

Now, writing down explicitly the Riemannian covariant derivative, we recast ()24p into 



{} 



(23) 



(24) 



(25) 



(26) 



Here, the first term on the rhs contains the full Riemann-Cartan connection, By - ! y 



K i:j k , cf. with EH 



It is also possible to write the conservation law in a different form. By raising the index i, we then straightforwardly 
can recast (l24l) into 



klj 



(27) 



Thus, equation (42) of [48[ is correct, it coincides with (|27|) . However, one should be careful since the position of 
indices in the definitions of the connection, torsion, contortion, and curvature is different from our conventions. Note 
also that the spin in Yasskin and Stoeger is defined with the | factor, see their definition (8) in |48j] , and compare 
it with our definition (IA4|) . It is satisfactory to see that our computations regarding the conservation laws are in 
complete agreement with those of Yasskin and Stoeger in [48| . 
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t = const 



/T'J 
t» j= 



/ A klj q\ 



FIG. 1: Sketch of the hypersurface E, i.e., the world tube of the test particle. A continuous curve through the tube is 
parametrized by Y a . Coordinates within the world tube with respect to a coordinate system centered on Y a are labeled by x a . 



IV. PROPAGATION EQUATIONS 

Let us switch to a notation which is close to the one in [48| . It turns out that (fT7|) is more appropriate to bring the 
energy-momentum conservation equation into a form analogous to the result (42) in [48]. By raising one index and 
explicitly rewriting 5 the covariant derivative in the first term of (jTTJ) , we obtain 

'/''•',, = />'.,/ /A"-' - Tkj'TW + N'Jv.A^ . (28) 
Furthermore, the hypermomentum conservation equation in (|2ip takes the form 

A klj d = N mj k A jlm - r m /A roW - T m j l ^ Hmj) - N jlm A k mj - f lk + V*. (29) 

By using (|21D in (117|) , we can also obtain an alternative version of (|28p , which has a very similar structure compared 
to (42) in {48|: 

l = R l 3k iA k[ 3 + N l kl N a k A^ la - N i kl N^ la A k aj - N l kl f lk - r fc /f + N\^ 1 

&T ij tj = RtjkiAW - N l H f lk - Tk/T^ + N i ki t kl . (30) 

Note that in the last equation R z j k i represents the full curvature. The structure of equation (f30)) is very similar to 
(42) in (48j . In the following we are going to derive the propagation equations for the integrated moments following 
from the conservation equations (|2"9"|) and ([317)) . 



5 Remember Vj (jS** + A'^ = -7= [\ r ~9 (S i: > + A*^)] . + V kj i S kj , where 5 ij denotes the symmetric and A^ denotes the antisym- 
metric part of a quantity with two indices. 
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A. Lemma: Derivative of the integrated moments 

The following relation, cf. (41) in [48| , between the time derivative of the multipole expansion of a current also 
holds within metric-afhne gravity 



(31) 



Here J a denotes the density of a matter current, in our case A kl i, T ij , or t kt . Additionally, we have Sx a :=x a -Y a , 
and p b a — 5x b , a = b b a — v h 5® = 5 h a — SqS° = 5^5" for the spatial projector. The upper-index of Ja u is associated with 
the last index of the corresponding matter current, e.g., Ja° —> T l °. In (f3"Tj) . and in the following, integrals are taken 
over a 3-dimensional slice E(f), at a time t, of the world tube of a test body. We use the condensed notation 



/= / f(x)d*x. 

E(t) 



B. Conservation equations integrated 

With the help of (|3"Tj) , we derive the integrated version of the conservation equations ([3(1 



d 
dt 



ib f3 _ y b f3 



J \a = l J 0=1 y J y a =l,a#/3 / 

f[ Sx b A (l<',,,A>' ! - N' l kl f lk - r ki l fw + x' u i u 



J] 5x h " ) T i0 



\a = l 



and 129 



d_ 

dt 



J \a=l / (3=1 \_ J y a =l,a^/3 / 



n Sx b ° ) A 

a=l,a=£j3 



A70 



With the introduction of new names for the integrated moments 



■bi—bnij 



the integrated conservation laws take the following form: 6 

v bf>T 



dj=bi—b n iO _ y-> f-=f>i-l -b n ib ^^bi—lp-bniO 



dt 



(3=1 



(32) 



Note that we use an inverted circumflex, e.g., bp, to indicate that an index is omitted from a list. 
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+ J (fl Sxbc ^j ^mj*A J,m - r ro j fc A my - r m /A fc (" 1 ^ - N jlm A k mj - f (fe + t 



(33) 



(34) 



Equations ([33]) and ([34"]) may be compared to (51) and (52) in (48[. 

C. Propagation equations for pole-dipole particles 

Let us now proceed along the lines of [48[ and derive the propagation equations for pole-dipole particles by using 

and (|3"4"|) . Here we investigate the case in which the following moments are nonvanishing: A ufc ,T u , , and t* 5 '* 
- i.e., we only take into account a pole contribution from the hypermomentum; the canonical energy- momentum and 
symmetric energy-momentum are considered to contribute at the pole as well as at the dipole level. This assumption 
is in accordance with the treatment in [48| . in which only pole contributions of the spin current were considered. Let 
us expand the geometrical quantities around the worldline Y(t) of the test particle, cf. figure [TJ into a power series in 
Sx a = x a - Y a . We have 

jkl\ x — R l jkl \ Y + fix" 1 R l jkl,a \ Y + ' ' ' J 
{} 



r ij 



r ^ 



N\i\ x = N l kl \ Y + Sx a N l kl , a \ y + ---. (35) 

The general form of the integrated conservation laws (l33|) and (|34|) then yields the following set of propagation 
equations: 

^T^ m ~\ U] ATi W lk ATi 7 T alk Q i7p( k i) U i 777 a ( fe i) 

— 1 = K jfejA — iV kl-L ~ kl,a-L — 1 kj 1 — 1 kj ,a-l 

+N i u t l +N i kl J akl , (36) 

jT l ° = T a - v a T° N* kl T lk - T k iT am + N i u r k \ (37) 

= T Ua + T alb - v a T m - v b T aM , (38) 

| A m = N mj k A jlm - T mj k A ml3 - T mj l A Hmj) - N^A k mj - T lk + 1\ (39) 

= A* ,o - t ,-A* ,0 -T° ,fc +r W . (40) 

Here we suppressed the dependencies on the points at which certain quantities are evaluated. The set (|3l)j) - (|4TJ|) 
represents the generalization of the propagation equations (63)- (67) in [48j to metric-affine gravity. 

V. ALTERNATIVE FORM OF THE PROPAGATION EQUATIONS 

It was pointed out by several authors, see also page 2086 in [HI, that the form of the propagation equations depends 
on the definition of the integrated moments, in particular, the index position in the set of equations (|32p . Of course 
ambiguities emerge due to the integration process and the fact that the metric is not a constant. In the previous 
section we used the index positions which match the ones used in [481 ] : this allows for a direct comparison of their 
propagation equations with our result in metric-affine gravity. Since there is a priori no way to tell which index 
position in the integrated moments is the more physical one, we are also going to derive an alternative version of the 
propagation equations, in which integrated moments with mixed indices are used. 

From a formal standpoint, the definition with mixed indices may be favored over the definition with upper indices. 
Geometrically, the momentum should always be a covector, i.e., it should have a lower- index. This becomes imme- 
diately clear if we recall some basic facts from classical mechanics. The velocity is a vector (with an upper- index), 
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v a = q a . Then, the momentum is, by definition, p a := dL/dv a - which obviously is a covector. Hence, from this 
standpoint it appears plausible to consider the choice 

p - It 

as definition for the momentum. In the following, we are going to work out an alternative set of propagation equations, 
which are based on the definitions with mixed indices. 

Once again, we start by rewriting the conservation equations for the canonical energy-momentum current (|17j) and 
hypermomentum current (|21[) . which take the following form 

TV j = iV^V + rVTy +N l3 k ¥ k , (41) 
AYj = r,z m A V - r m /AV l -f , fc +?i. (42) 

Note that Tij k represents the full connection, the last two equations should be compared to (42) and (43) in [48| . 
Apart from the index positions, equations (|41 [) and (|42f are completely equivalent to (|30j) and ([29]) . 

A. Conservation equations integrated 

With the help of IpTl]) , we derive the integrated version of the conservation equations (I4T]) : 



and ([291) 




Now we introduce the integrated moments with mixed index positions. Note that we use an underline (lower-index 
position) to distinguish these definitions from the overlined (upper-index position) quantities in (|32p 




With these definitions the integrated conservation laws take the following form 




11 



J 



d 

A" 1 i" = 

13=1 

+ j(fl 5xba J (V"^ V ~ T m /A^ m - f t k + ?,) . (45) 
Equations (|44ll and (j45| should be compared to (|33|) and (f34|) . as well as to equations (51) and (52) in (48|. 

B. Propagation equations for pole-dipole particles 

Finally, we derive the propagation equations for pole-dipole particles by using (|44|) and (I45j) . Again we investigate 
the case in which the following moments are nonvanishing: A 1 j k , Tj , T 1 j k ,t l j, and t y & . The expansion of geometrical 
quantities around the worldline Y(t) of the test particle, cf. figure HI into a power series in Sx a — x a — Y a , reads 



Rijk 


a; 




i 1 


-j- (5x Rijk ,a | y 


"P & 
1 V 1 


a; 


P & 1 
1 *j | 




(5x r.jj ' ,a | y H" 


AT- 


X 


Nii" 


\y " 


K<5z Q % fe ,a| r - 



(46) 

The general form of the integrated conservation laws (|4~4"|) and (|4"5j) then yields the following set of propagation 
equations: 

j t Z° = iWAV + T^Tj + r,/, a r fe J + NijHh + N l3 \ a t a h, (47) 
d_ z aO = T i a -v a T i °+r ij k r t k i + N fj k ?* k , (48) 
= T b l a + T a . b -v a T b l °-v b T a i °, (49) 



-^aY = iy n A V - r mj k ^r - r t k + t k h (50) 



= A?-/AV-TY+f|. (51) 

Again we suppressed the dependencies on the points at which certain quantities are evaluated. The set (|47[) - (|5T]) 
represents the generalization of the propagation equations (63)-(67) in [48j to metric-affine gravity, now with the 
mixed index convention. The above set should be compared to our result in (f56")) - ([4T))) . 

C. Rewriting the propagation equations a la Yasskin and Stoeger 

Now let us rewrite the propagation equations of metric-affine gravity (|47|) - (|5ip in a form which closely resembles 
the main theorem of Yasskin and Stoeger in Poincare gauge theory, i.e., equations (53)-(58) in (48|. We start with 
the following identity which holds because of the definition of the projector //V 

Ab a a \b i a \b k /rn\ 

_ c = V A c + p k A c . (52) 

Using this relation, the last one of the propagation equations (f5Tj) takes the form 

T alk -t akl = p a b A klb . (53) 
This equation may be compared to equation (68) in [48l] . Again with the help of (|5"2"j) we can rewrite ([50)) as follows: 

t k t - T k = V„A fc j° + (r jm k A m i b - r/A 1 ,, 1 ) ph, (54) 

where 



V„AY := ^A fc ,° + v "T m /AV - v m T ml 'A k ] °. (55) 






rpbia 


|_ rpaib 


— y a T bi ° 


— y b T ai ° 





rpiab 


■ rpbai 









rpabi 


■ rpiba 


_ y irpabO 


_ y arpibO 
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Equation (|54[) should be compared to equation (69) in [48[. Proceeding along similar lines as in we are now 

going to cyclically permute the indices in (|49p twice, resulting in 

(56) 
(57) 
(58) 

Then adding ([56")) and ([57)) and subtracting yields 

Q _ rpb(ai) _ rpa[bi] _ jri[ba] _ y arp[bi]0 _ v brp(ai)0 _ y irji[ba]Q /ggs 

This equation should be compared to (72) in [43]. We proceed with the following identity: 

^rp(ab)0 22" lQ [^] _j_ 22^[ a 0] -j- _J_ r p^ a (60) 

Combining (|60p with (|49|) . in which we raise the index and set i = 0, we arrive at 

2ji(ab)0 _ 2y a [W] _|_ 2 r p b l aQ ] _|_ ^ay^OO _|_ ^fcyaOO (61) 

which should be compared to (74) in [48|. Remembering that 

2rp[a0]0 _ rpaOO (62) 

which follows directly from Sx° = 0, we can rewrite (|61[) as follows: 

y(afc)0 = v (a^b)0 + p a m ^[0b]m + ^[06] + ^^[0"]". + ^[0a] > (63) 

where we made use of (|53j) and introduced the following definition for the antisymmetric part of the integrated orbital 
momentum on the basis of the canonical momentum 7 

A ab := 2T [abI °. (64) 

Remembering that t ab is a symmetric quantity, equation (|63p can be rewritten as 

T (ab)0 = i/ a A b >° + p a m Aj ob l m + p b m A[° a ] ,n , (65) 

which is analogous to equation (76) in [4^]. This result can be used to rewrite ([59]) . 

r^b(ai) rpa[bi] _j_ rpi[ba] _j_ ^arj->[bi]0 _j_ ^irp[ba]0 

which resembles the first part of (77) in ^48| and can finally be brought into the form 

T b{ai} = p b m [-v {a k i)m + p a n A [tm]n + pVAl™'") , (67) 
which is analogous to the second part of (77) in j48|. The last equation can be used in (|53[) to obtain 

t akl = p a b (A (fei)b + v^ l k k)b - p l n A^ n - p k n A^ n ^j . (68) 
After reinsertion into ([53)) we arrive at the final result, 



T alk = p a b fl^kb + ^kib + v (i^k)b _ p l n ^[kb]n _ p k n ^lib]n\ ^ 



(69) 



7 Note that this definition corresponds to the quantity L ah in |4q| . In this work, in contrast to [48| . we use the symbol L ab for the 
"complete" first moment of the integrated canonical momentum, i.e., including also the symmetric part. 
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which closely resembles the form of one of the propagation equations found [48j], i.e., equation (56). With the help of 
(fSU)) . (f6"5|) , and (f6"9")) . equation (|48[) can now be transformed into 



d 
~dt 



+ 9u (V a A' )0 + P a m A^ m + p l m A l ° a] 
lyfc A Qa/^ + A jkb + v^A^ - p k n A^ n - p> n A} kb ^ 
N ijk p a b (£S 3k)b + i/ fe A j)fc - p k n A [3b]n - (P n ^ kh]n ) , 



(70) 



where we introduced := T,- for the integrated 4-momentum. Equation (|70[) is analogous to the propagation 
equation (55) in With the help of (f70|) we can can bring (f54|) into the form 



N ljcP k b (A« c ) b + t/ c A^ b - p c „Ab' b l" - p J „A[ cfc ]") , 



(71) 



which can be viewed as the analogue to (79) in [48j| . Because of the different symmetries in metric-affine gravity 
the method used in this section, which was outlined in (48|, does not lead to a very compact form of equation (54). 
The last equation in the rewritten set is the one relating the time derivative of the momentum to the other matter 
quantities; from (gTJ), (O, and (JHSJ) we obtain 



d 
df 



Ei = R i jk l (v j A k l + ft n A k S 



+V [v„aV + p n t (r n J A m k l - T nk m A\ 

~t~r^" \aP k r ij i± k H" r ij .at. ^ k- 



(72) 



This equation should be compared to (80) in [48|]. We only note that an elimination of P k and t a ^ k in the last two 
terms of (|T2"|) is possible by using (j53"|) and ([70]) . In the next section we work with a slightly different set of quantities, 
which allow for a very condensed form of the propagation equations of metric-affine gravity. 



D. Rewriting the propagation equations 



In this section we present a more condensed form of the propagation equations. Thereby we find a direct general- 
ization of the main result 8 of (48|, i.e., equations (53)-(58), in the case of metric-affine gravity. 

We introduce the following notation for the integrated quantities: :— T_® denotes again the integrated 4- 
momentum, L fe ; := T k i° the total orbital canonical energy-momentum and Y_ k i := A fc ;° the integrated intrinsic 
hypermomentum. Furthermore, recalling that the hypermomentum comprises the spin, dilaton charge, and intrinsic 
shear, it is convenient to denote the antisymmetric part of the hypermomentum as the integrated spin T k i :— A^ fe ;]°, 

whereas the trace of the hypermomentum defines the integrated dilaton charge Z_ := A k k °. 

In addition, we introduce a shorter notation for the "convective currents," i.e., the projected quantities which we 
have used in previous sections and which are also used in (48] ] . For the intrinsic hypermomentum, we have 

(c) 

Ak m . \km „,m A k — „m \k n 

and for the orbital canonical energy-momentum 

(c) 

Tk m rrik m „,m rpk rn rpk n 

_ I ■= i_ 2 — V J_ I — p n ±_ I . 



Please note the typo in equation (53) of |48ll . Using the notation of [48ll the last term in (53) should read: . . . + vN^ a,y g Je S/ e \ a gs. 
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The convective spin and dilaton currents arise as the antisymmetric part and the trace of the convective current of 
the intrinsic hypermomentum, i.e., as 

(fl'ft m . a \k m „,m k 

z i ■= a i] - V T I 

and 

(c) 

Z_ k := Z k -v k Z, 

respectively (here Z fe '■= A-*j k ). With this notation, we recast the propagation equations (|4"5|) - ([51"j) into 



Z k l = v l P k + j t L\ - Tk/Ti 3 + N k / A V, (73) 

(c) 

T^^ = 0, (74) 

V„r fc = —T_ k l +t i k - Tjf AV + Tjk 1 AV, (75) 

A fe / a = T a i k -t ak i. (76) 



Equation (|73[) describes the canonical energy-momentum in terms of the usual combination of the "translational" 
plus "orbital" contributions (the first two terms), plus the additional contribution of the first moments. One should 

(c) 

compare this with the alternative formula (|7D|) . Equation ([74]) simply tells us that the convective current T_ a i is 
antisymmetric in the upper indices a and b. This is a useful technical fact. The next equation (|T5[) is actually an 
equation of motion for the intrinsic hypermomentum. Its form closely follows the Noether conservation law of the 
hypermomentum, cf. (|19p and (|21[) . An alternative form of such a dynamical equation for the hypermomentum is 
given in (|71[) . Finally, the equation (|76[) expresses the convective intrinsic hypermomentum current in terms of the 
first moments of the energy-momentum. 

Equations (f73")) -(|76 p are easily derived from (j4"51) - (f5T]) . one only needs to rearrange some terms. In contrast to this, 
we need some additional steps to arrive at a new form of equation (|47p , which represents the most interesting of the 
propagation equations from a physical point of view. 

We start by expanding the general connection in (|4"T|) , this yields 

^ = iWAV + { lk%* - N tk > (T k - tS) + Tik l ,aT a i k ~ N lk \ a (TV - ?\) ■ (77) 
Furthermore, we have 

j t (l? - A^A fe ; °) = ^T, - v a N ik \ a A k t° - N lk l j t A k i°- (78) 
Insertion of ([50]) and dZTj into flTSJ) yields 

-r mj k A j r) - N lk l ,„ (T /" - t ak i + v a A k l °) (79) 

= Tik% k + Tik l , a X_ a i k + A k S {R ljk l - T 3P l N lk P + T jk PN ip l 

-N ik l d + T Tl p N pk l - T^N pk l \ . (80) 

In the last step we made use of (|5T|) in order to replace the terms in the last brace in the second line of (|79p . 
Furthermore, we added a "0" dummy term, i.e., the last two terms in the second line of ([8H)l . We proceed by replacing 
the curvature by its decomposition, i.e., 

Ri jk l = Rijk 1 + VjNi k l - ViN jk l + N v l N ]k p - N jp l N ik p , (81) 
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equation (180)) then turns into 

d /„„ „ , . ; ,A U - U , _.: , .,;/!! , U , Li 



| (T, - ATifc'A*, ) = r ik % k + Tik l ,j Zh k + A V ( R ijk l ~ Tji p N pk l - V, V ) . (82) 



We rewrite (j48J) with the help of (j51j) : 

= jT. k i° + « fe Tj° - Ti P m T k m P + N lp m (A p m k - v k A p m °) . (83) 
Contracting this equation with the Levi-Civita connection and introducing another "0" dummy term yields 



i i T"^ _ ^ m T 1 p _r_ / AT m a p k 

1 ifc ,a±_ / L ifa t lp ±_ rn T" 1 ik -*-*lp ±± m 

rV rV n - rV rik m ) U> m v k . (84) 



With this result at hand we can replace the first term on the rhs of (f82|) . i.e., 

= (r.i.'j - ",/",,'") i V + aV (r«*' - v.iV) + Rm'I'.V. (85) 



If we introduce a new quantity 



Vi := T° - N ik l A k l ° - T lk l T k l ° (86) 



as a generalized total 4-momentum, equation (|85p can be written in a more compact form as follows: 

v„7>< = ^Rkji l T k i° v j + (r ik l j - Ti* rV) 2 j i k + AY (izW - vW') . (87) 
{},{},{}, 

By using the Ricci identity Rjki + Rkij + Rijk = and the fact that the convective part of first integrated moment 
of the canonical-momentum is antisymmetric in the upper two indices, i.e., T m p =T_ ' m p ', we can recast (|87[) into 

V v Pi = R kj i l T k l °v j + ( { R ijk l - ViNjA AY + { R ljk l T V- 



This equation represents the rewritten form of (|47j) and should be compared to (|72|) in the previous section. 

It is worthwhile to notice the general feature that characterizes the coupling between the physical objects (currents) 
with the geometrical objects (metric, connection, and the derived quantities). Namely, the intrinsic current (the one 
that is truly microscopic, which arises from the averaging over the medium with the elements with microstructure, i.e., 
that possess internal degrees of freedom) couples to the post-Riemannian geometric quantities, see the second term 
on the rhs of ([56")) and the first term on the rhs of ([55)1 . In contrast to this, the orbital canonical energy- momentum 
(which is induced by the macroscopic dynamics of the rotating and deformable body) is only coupled to the purely 
Riemannian geometric variables and never couples to the post-Riemannian geometry, see the last terms on the right- 
hand sides of (|86|) and (|88| . This observation represents a generalization of the result of Yasskin and Stoeger [48| , in 
other words, it proves that the possible presence of the post-Riemannian geometry (in particular, of the torsion and 
the nonmetricity) can only be tested with the help of the bodies that are constructed from media with microstructure 
(spin, dilaton charge, and intrinsic shear). Test particles composed from usual matter, i.e., without microstructure, 
are not affected by the post-Riemannian geometry, and they thus cannot be used for the detection of the torsion and 
the nonmetricity. 

In order to get a better understanding of this fact, we will consider several special cases of the metric-affine geometry 
in the subsequent sections, moving from a general non-Riemannian geometry back to the Riemannian one. 
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VI. RELATION BETWEEN THE INTEGRATED MOMENTS 



In different situations, it is technically convenient to use different definitions of the integrated moments (see also 
[85j for the behaviour under infinitesimal coordinate transformations). However, directly from the definitions (|32|) 
and (|43p we can establish relations between two sets of the moments. 

Starting with the identity = g 3 we expand the metric in the same way as the other geometric quantities 
@6J>, 

9 jk \ x = g ]k \ Y + 6x a gi\ a \ Y + ■ ■ ■ , (89) 
and then by integration over the world tube, in the pole-dipole approximation we find 

t* = f j - 2 Ti {kj) tk- (90) 

We used here the metricity condition g^ k ,i = — V h-J g nk — T i n k g jn - 
Analogously, we have for the integrated canonical energy-momentum 

T ij =T ij -2rV ,fc) TV. (91) 

The "inverse" formulas read 

tj = r j+ 2r lUk} t Hk , (92) 

Xj = ^+2r ((lt) f J . (93) 

Hence, in the pole-dipole approximation, the integrated hypermomenta and the first moments of the canonical and 
metrical energy-momenta in both sets are the same: 

A ijk = A l3k , (94) 

ty k = t jk , (95) 

T ijk = T ]k . (96) 



With the help of |95|) and (|96|) . we can verify the consistency of the relations (|90|) and (|92|) . as well as (|9ip and (|93p . 

For single-pole test particles, the corresponding integrated energy-momenta coincide since the last terms in (|90l) - (f9"3")l 
vanish. 



VII. CONCLUSIONS & OUTLOOK 



In this work we derived the equations of motion for test particles in metric-affine gravity from the conservation 
laws of the theory with the help of a multipolc formalism. Apart from the general form of the equations of motion, 
we explicitly presented the propagation equations for pole-dipole test particles. Our results are valid for a very large 
class of gravitational theories, i.e., all theories which fit into the framework of metric-affine gravity. The equations 
derived in this work should be used to systematically study the motion of test particles with spin, shear, dilation, 
and rotation within alternative gravitational theories in a non-Ricmannian context. Our results generalize previous 
analyses [HI [IB, [13, [H, which were carried out in the context of general relativity, Einstein-Cartan theory, and 
within Poincare gauge theory 



A. Special cases 



In this section we discuss several special cases within our framework by either making assumptions about the internal 
structure of the test particles, or by constraining the background geometry. The full agreement, in some special cases, 
with the well-known results from general relativity and Poincare gauge theory demonstrates the consistency of our 
framework. 
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1. Equations for a single-pole particle in metric-affine gravity 

Let us consider the propagation equations for a single-pole test particle in metric-affine gravity, i.e., the set (j36j) - (j40)l 
with vanishing dipole contributions: 

v a T l ° = T*°, (98) 

| A H ° = AT ro /A'' im - t/^ - r m /A k(m3) - N^A k mj - + 1\ (99) 

v a A kW = A kla . (100) 

It is a common folklore that in generalized gravity theories the equation of motion for single-pole test particles is given 
by some kind of "generalized" geodesic equation. By generalized we mean an equation which has the same form as 
the geodesic equation, i.e., the equation of motion for single-pole test particles in general relativity, but in which the 
Levi-Civita connection has been replaced by the full (non-Riemannian) connection. The result in (f97|) - (ll00j) clearly 
demonstrates that such an assumption is not substantiated. 

a. Particles without intrinsic hypermomentum If we perform a further specialization by considering only test 
particles without intrinsic hypermomentum, the set (|9T[) - (jl00[) turns into 

|r° = -tfvr* - { A/t (w) + n*j*, (ioi) 

v a T° = T a , (102) 

T lk = t l . (103) 

Of course the first and the last term on the rhs of (|101[) cancel because of (|103[) and the equation of motion for a 
test particle without intrinsic hypermomentum is then given by the regular geodesic equation [in the next section we 
explicitly show how one can recover the geodesic equation from the set (|101j) - (|103[) ]. This generalizes the well-known 
result from Poincare gauge theory to metric-affine gravity, i.e., a test particle without intrinsic hypermomentum will 
not "feel" the torsion or the nonmetricity of the underlying spacetime. Hence, test particles without intrinsic spin, 
shear, or dilation current are not suitable for mapping the non-Riemannian features of spacetime. Accordingly, current 
experiments like Gravity Probe-B [9(| are not suitable for the detection of torsion in contrast to what is sometimes 
claimed by other authors. At this point, one should mention that a coupling between torsion and matter without 
intrinsic spin currents may be achieved in some nonstandard gravity theory, although the authors of the present 
paper are not aware of any viable candidate for such a theory. For any theory which fits into the very general and 
well-motivated framework of metric-affine gravity, e.g., Poincare gauge theory and Einstcin-Cartan theory, such a 
coupling will not occur. 

2. Recovering the geodesic equation 

In this section we explicitly show that the single-pole equations of motion for a test particle without intrinsic 
hypermomentum take the form of the usual geodesic equation. The set (I10ip - (|103[) reduces to 

jT W = - { r } fe /T W) , (104) 

T a = v a T , (105) 
T lk = t\ (106) 

Now lets us introduce the velocity u a := dY a /ds along the world line of the particle. Note that u° — dt/ds, 
ds 2 = g ab dY a dY b , and remember that Y a (t) = x a (Y(t)) = t6$, d/dt = v a d a , u a u a = 1, v a = dY a /dt. With this 
definition we can rewrite (|104|) and (|105[) as follows: 

ArVr fe >°2* = 0, (107) 
ds 

u°T la = u a T l °. (108) 



18 



Setting i = in the last equation and reinsertion into (|107j) . together with the definition m := T°°/ (u ) 2 , yields 
T* a = mu l u a . This in turn can be used to rewrite (|107[) as follows: 

d_ 

ds 



(mu l ) + T kj i mu k u> = 0. (109) 



Multiplication of this equation by u.- L and remembering that u h V&u° — (u a ^ + T c b a ' ! ' j du"/ds (/".;,//'' 
« a Vf)M" = yields 

dm ■ . O k dm 

——uUi + mukuS/iU =() => — — 

ds ds 

When we use this result in (|109[) we end up with 



^-rfin + mu k u j VjU k = => -j- = 0- (HO) 



^- + r«W - o, (in) 

ds 

which is the geodesic equations. Hence, in metric- affine gravity single-pole test particles without intrinsic hypermo- 
mentum, i.e., without spin, shear, and dilation currents, move in exactly the same way as test particles in general 
relativity. We stress that no constraining assumptions about the geometry of the background spacetime have been 
made in order to derive this result. Equation (jllip is valid in a completely general metric-affine spacetime, i.e., the 
background can be a non-Riemannian one with nonvanishing torsion and nonmetricity, the test particle just does not 
feel this geometric features as long as it does not posses any "microstructure" in the form of a nonvanishing intrinsic 
hypermomentum. 

In later sections we are also going to discuss the equations of motion for some special cases in which we impose an 
a priori restriction on the geometry of the background spacetime. 



3. Recovering the Mathisson-Papapetrou equations 

Also the well-known propagation equations for a classical pole-dipole test particle can be easily recovered in our 
framework. For particles without intrinsic hypermomentum in a Riemannian background the propagation equations 
in ([36 |1 -(|40 ]) turn into 

p M = ^-If^, (112) 

d_ T aiO ^ T a_ vaT _ U k . ir m (n3) 

v-T m +v b T ai ° = T bla +T alb , (114) 
¥ k = t\ (115) 
T lk = t akl . (116) 

These equations are exactly the equations of motion for a pole-dipole particle described by Papapetrou in (3.2)-(3.4) 
of [2l|. This result clearly demonstrates the consistency and generality of our framework. 



4- Propagation equations in a Weyl-Cartan spacetime 

The Weyl-Cartan spacetime is characterized by a special type of nonmetricity, namely, when the 1-form of the 
nonmetricity Q a p — g a p Q reduces to just the Weyl covector Q = Qidx 1 . Correspondingly, the distorsion 1-form then 
reduces to 

Nj^-^iQ + Kj, (117) 

where the contortion K a p = — Kp a := N\ a m is just the antisymmetric piece of the distorsion (note, however, that 
KoP is constructed from both the torsion and the Weyl nonmetricity). In components, we have explicitly NiJ 3 = 
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Using relation (|117[) , we derive the propagation equations for test particles on the background of the Weyl-Cartan 
spacetime: 



{} {} 



R kj i l T k i°v j + ^R ijk 1 - Vi-Ki*') zV + Rijk 1 TV + \$iQo)Z j , (118) 

= v l p k + j t L\~ rV ry + K kj l - \Qk k\ (H9) 

(c) 

T (o i 6 ) = 0, (120) 

V» r fc = - ZV + f fc - r^i* 4V + T jk l A'V + A.V - K- jfc J AY, (121) 

(c) 

AY = T a i k -t ak i. (122) 

Here V i =P i + \Q i Z- K ik l I k l - T ik l L k h 

a. Single-pole particles For the single-pole case (when all of the first integrated moments vanish), we find a 
surprisingly simple system 



VvP z + K tl k viP k = R ijk l v jT k l + -f ijV iZ--Q i ^, (123) 



£} „ „ , ,„ „ , , u 1 j 1 dZ 

- JijV J Z — -Qi —r- 

2 J 3 — 2 dt 
T k = v l P k , (124) 
V V Y\ = -T^+tk, (125) 

A k i a = 0. (126) 

Here we introduced /y := diQj — djQi. Thus, provided a test particle has a nontrivial integrated dilaton charge Z_, it 
will be affected in the Weyl-Cartan spacetime by the Lorentz-type force represented by the second term on the rhs of 
the propagation equation (| 1 23[) . If, in addition, the test particle has a nontrivial spin r k i, the latter will be affected 
by the Mathisson-Papapatrou-type force which is determined by the Weyl-Cartan curvature, as described by the first 
term on the rhs (|123[) . 



5. Propagation equations in a Weyl spacetime 

Weyl 0, |aE SI was the first who noticed a similarity between the electromagnetic vector potential and the non- 
metricity covector Qi. Indeed, this is also manifested in the equations of motion, as becomes clear from the rhs of 
equation (|123p . However, an essential difference is that the Weyl nonmetricity may interact with the dilaton charge 
and not with the electromagnetic charge. 

The Weyl geometry arises as a special case of the Weyl-Cartan spacetime, when the torsion Sij k := Tij k — Tj k — 
is equal zero. In this case the distorsion is still given by (|117|) . but the contortion is expressed in terms of the Weyl 
covector only: 

K ij k = ±(g ij Q k -6>?Q j ). (127) 

The propagation equations in the Weyl spacetime are formally the same as (|118l) - (|122[) where we have to substitute 
the contortion (|127ll . Analogously, the dynamics of single-pole test particles is described in the Weyl spacetime by 
(P3 > (fT251 ) with (TT27)) inserted. 



6. Propagation equations in a Riemann-Cartan spacetime 

The Riemann-Cartan spacetime arises from the Weyl-Cartan geometry for the case of vanishing nonmetricity, 
Qi = 0. The distorsion then coincides with the contortion and is constructed only from the torsion: Nij k = Kij k — 

2 {Sjki "I" Si k j + Sji k ). 



9 Our notation for the spacetime torsion is different from [5dl |. Since we reserved the symbol T for energy-momentum related objects, the 
torsion tensor is here denoted by the symbol S as in the old review [9^1 . 
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= u\iP jV 




= u ° (%) " 



The propagation equations for pole-dipole particles in Riemann-Cartan spacetime are easily derived by putting 
Qi = in equations (|118[) - (|122[) . We will not write these equations explicitly. 

a. Single-pole particles In order to discuss the propagation equations for single-pole particles, we again introduce 
the 4-velocity u a :— dY a /ds along the world line of the particle. With u° — dt/ds and ds 2 — g a bdY a dY b , we have 
u a u a = 1 (note that u a — u°v a ). Then, it is straightforward to verify that in the Ricmann-Cartan spacetime equations 
(fl23 )) -(p5 | reduce to 

±Ri jk l u j T k h (128) 

(129) 
(130) 

Tm) • ( 131 ) 

Here we denoted the covariant (Ricmann-Cartan) derivative along the trajectory by a dot: "'"= D/ds = u l Vi. 

It is satisfactory to see that with (|128[) and (|130p we recover the usual equations of motion for a test particle 
with mass and spin in the Riemann-Cartan spacetime [86l [87], |93| . One should note, however, that we are still 
in the framework of the metric-afhne gravity in which a test particle carries, besides the mass and spin, also the 
dilaton charge and the intrinsic shear. The latter integrated characteristics are described by the symmetric part 
of the intrinsic hypermomentum Y_{ij) . The dynamics of these quantities is determined by equation (|131j) which is 
completely decoupled from the other propagation equations. In other words, they do not affect the motion of a test 
particle in the Ricmann-Cartan spacetime, and the trajectory is completely defined by the behavior of the integrated 
4- momentum and the integrated spin r k i. 

Let us contract equation (|130[) with u % . This then yields the explicit form of the integrated 4-momentum, 

Pj =muj +2u i ± ij , (132) 

where we introduced the notation for the rest mass of the body m := (i.e., the momentum projected to the rest 
frame). By substituting this back into (|130jl we obtain the dynamical equation for the spin 

±ij - UiU k ± k3 + UjU k ± ki = 0. (133) 

7. Propagation equations in a Riemannian spacetime 

When all the post-Riemannian geometric objects are trivial (no torsion and no nonmetricity, i.e., N a @ = 0), the 
propagation equations on the purely Riemannian spacetime reduce to 

(c) , 

(134) 
(135) 

(136) 
(137) 
(138) 
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a. Single-pole particles For the single-pole particles with vanishing intrinsic hypermomentum this simplifies to 

VvE l - 0, (139) 
ZC = v { P k , (140) 
T-k - tk- (141) 



The resulting trajectories are geodesies. 
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8. Propagation equations in a Riemannian spacetime (alternative form) 

For completeness let us also determine the explicit form of the propagation equations using the upper-index con- 
vention for the integrated moments, for the special case of a Riemannian background. From (|36p -(|40 p we can infer 
that pole-dipole particles move according to 

^7-10 O i -r-klj O .— (kj) ,„ 

-^-T — R \ ]U A - r kj T , (142) 

v a T l ° = T ia , (143) 

— A = —T m j A - T m j A J '-T +t , (144) 

- A Ma - v a A m - T lk + T kl . (145) 

a. Single-pole particles Further restriction to single-pole particles with vanishing intrinsic hypermomentum A abc 
brings the set (|142[) - 1|145[) into the form 

dt T = - Tk > T ' 

r a = v a r\ 
T lk = t l . 

As we have already shown these equations lead to the geodesic equation. Hence, within our general formalism we can 
quickly reproduce the standard result of general relativity. 



B. Open problems 

The results obtained in this work are valid for a wide class of extended gravitational theories that are naturally 
embedded into the framework of metric-affine gravity. However, our study is not exhaustive in many important 
aspects, and at this stage there remain several interesting open questions related to the multipole expansion of the 
equations of motion of test particles in alternative gravity theories. 



1. Invariant definition of moments 



As we have already mentioned in previous sections, the definition of the integrated moments of the matter currents 
in the multipole formalism is to a certain extent ambiguous. This is related to the index positions in the integrand 
expression and to the nonconstancy of the metric which is used to lower and raise the indices. In view of this problem, 
we decided to present the full set of propagation equations for two different choices of the integrated moments, defined 
in equation ((32]) and (|43|) , respectively. Thereby one covers the definitions which have been discussed most frequently in 
the literature. Although we clearly favor the definition with mixed indices (I43p . for the formal reasons given in section 
IVl even other index positions than the ones investigated in the present work are imaginable. Such an ambiguity 
in the definition of the integrated moments motivates the search for an invariant formulation. The corresponding 
program was already carried out in several works within a general relativistic context [3l|, [H, HH, HO, [94j • Within 
an alternative gravity theory like metric-affine gravity, which is no longer a purely metric theory but has a richer 
geometrical structure, a detailed investigation is needed in order to generalize the concepts linked to such an invariant 
formulation. 



2. Supplementary conditions 

Previous analyses [U [95|, [9(| [9?| in metric theories of gravitation have shown, that even at the dipole level, 
supplementary conditions are needed in order to obtain a closed set of propagation equations. Indeed, let us recall 
the propagation equations in the Riemann-Cartan spacetime, for example. The four equations (| 132[) are sufficient 
to find the four coordinates of a position of a particle on its trajectory. However, the system (|133|) contains only 
three independent equations, and this is not sufficient to determine six components of the spin. As a result, the 
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supplementary conditions are usually imposed on the spin of the test particles in order to make number of the equations 
equal to the number of unknown variables. The imposition of an additional supplementary condition comes with some 
assumptions about the physical nature of the particles under consideration, and there is no unique prescription how 
to do it. Even within the context of general relativity, a number of competing conditions exists. Furthermore, there 
seems to be no consensus on which of the supplementary conditions is the most physical one. In the context of 
alternative gravity theories the spectrum of possible supplementary conditions is greatly enhanced. This fact can 
be ascribed to the additional degrees of freedom within such theories, in particular, regarding the matter variables 
describing the internal structure of particles. Although there exist several studies of such supplementary conditions in 
the literature, most of them in the context of Einstein-Cartan and Poincare gauge theory, a systematic and up-to-date 
analysis in the context of metric-afhne gravity is still an outstanding task. We only note that an ultimate judgment 
over the correct choice of a supplementary condition can only be made with the help of an experiment. 

3. Propagation equations involving higher moments 

If we take into account previous results in Einstein's theory [98j , it is to be expected that the role of supplementary 
conditions is even aggravated at higher orders of approximation. Of course this is due to the fact that at higher orders 
we need an even more detailed description of the internal dynamics of the test particles. Nevertheless, the study of 
higher orders of the propagation equations, beyond the pole-dipole level, will be of great interest in the context of 
radiation phenomena. In particular, we expect that such studies will shed light on our understanding of the new field 
strengths of metric-affine gravity, i.e., torsion and nonmetricity, which have no counterpart in the classical theory 
gravitation, namely general relativity. 

4- Relation to other approximation schemes 

From a more formal standpoint, we can also ask about the compatibility with other approximation schemes which 
were employed in the context of gravitational theories. The most prominent examples being the post-Minkowskian and 
post-Newtonian approximation. Since these approximation schemes, in their full generality, are still under construction 
in the context of metric-affine gravity, a systematic comparison with the results obtained within a multipole scheme 
appears to be a long term project. 

To sum up, the study of the propagation equations of deformable test particles with the help of a multipole 
approximation scheme is a very rich field of research. In the context of alternative gravity theories this field is still 
in its infancy. Apart from the first steps undertaken in this work a number of open problems remain; we intend to 
attack these in future works. 
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APPENDIX A: GENERAL CONVENTIONS AND NOTATIONS 

In the theory of metric-affine gravity, the gravitational field is described by the three basic variables: the metric 
g a p, the coframe , d a , and the linear connection Tj 3 . The Latin indices . . . are used for local holonomic spacetime 
coordinates and the Greek indices a, /3, . . . label (co)frame components. The vector basis dual to the frame 1-forms 
■& a is denoted by e a and they satisfy e Q J$ /3 = 5@. Here J denotes the interior product (contraction) of a vector with 
an exterior form. Using local coordinates x z , we have $ a = hfdx 1, and e a = h l a di. All objects and equations that 
carry the local Lorentz indices can be recast into their counterparts with the coordinate indices with the help of the 
contraction with the components of the tetrads, hf and h % a . 
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1. Geometrical objects 

The geometry of MAG is described by the curvature 2-form R a " := dT a @ + A L Q 7 , the nonmetricity 1-form 
Qaj3 := ~Dg a p, and the torsion 2-form T a :— Dfi 01 which are the gravitational field strengths for linear connection 
r a ^, metric g a p, and coframe 1?", respectively. 

It is convenient to define a 1-form tensor-valued difference of the Riemannian (Christoffel) connection and the 
general linear connection: 

N J:=tJ-tJ. (Ai) 

This quantity is known as distorsion 1-form. In particular, the torsion is recovered from it as T a = —Np a A I?' 3 , 
whereas the nonmetricity arises as Q a p = —2Nr a py The corresponding curvature 2-forms are related via 

Rj = { RJ - DNj + A N a \ (A2) 



2. Physical objects 

The sources of the metric-affine gravitational field are the 3-forms of the canonical energy-momentum and hyper- 
momentum. They are defined by the variational derivatives of the material Lagrangian 4-form L mat , respectively: 

S« = %?. (A3) 
A% = ^f. (A4) 
The Lagrangian L ma t also depends on some matter fields ip, but this is irrelevant for the current discussion. 



3. Components 

When the local coordinates x l are chosen, we can write all the geometrical and physical quantities explicitly in 
terms of their components: 

(A5) 
(A6) 
(A7) 

(A8) 

(A9) 
(A10) 

Here r\ is the volume 4-form. Writing the Lagrangian form as L mat = £ ma t dx° A dx 1 A dx 2 A dx 3 , we can recast the 
definitions (|A3j) and (|A4|) as follows: 

Tj = (All) 
AV = (A12) 





= hfdx 1 , 


r 13 


— r @ dr 1 


N 

1 v a 


-N~ia. dx ; 


R 13 


— ~ q. dx f\ dx ^ 




= Tjdifr, 


A^ 


= AV^h- 



APPENDIX B: DIMENSIONS & SYMBOLS 



In order to fix our notation, we provide some tables with definitions in this appendix. The dimensions of the 
different quantities appearing throughout the work are displayed in table [TTJ Table IIIII contains a list with symbols 
used throughout the text. 
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TABLE II: Dimensions of the quantities within this work. 
Dimension (SI) Symbol 



Geometrical quantities 


1 


gap, 5 a p, g tj , ^J~^g, hf, T a f3 


, Nj, Kj, Rj, Q aP , Q, fp 


m 


x\ dx\ ds, Sx\ Y a , iT, T a 




m" 1 


P V- & N- P K ® O a 
Cot, L IQ 5 Juq 5 -ivja , l^mp 


Qi, Sij 


m- 2 


D.. f.. 
J^ija , Jij 




m 4 


V 




Matter quantities 


1 


u a , v a , p\, V 




kgm 2 /s 


h (Planck constant), L, L mat , L tot , A a fj , T a/3 , a " 3 , m " 3 , M af3 , H a p , E a p, A 1 /, A yfc , T l j k , t ij k , 




T l3k ,T k ,Y\,T\,L k u k* 


, r^i, Z k , Z 


kgm/s 


TT Tp V* rp k ,i rp % 3 T 

11 a , J-Ja. i Z-i a , J_| , b_ j , _f , L 




kg/(m s) 


A 




kg/(m 2 s) 


T a , C mat 




Operators 


1 


d, D 






{} 




m- 1 


di, Vi, V„, L 5 
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TABLE III: Directory of symbols. 



Symbol 




Explanation 


Form degree 


Differential form 


Component 






Geometrical quantities 




gab 


Metric 
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Determinant of the metric 


u 


V 




Volume torm 
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■& a 




Coframe 
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Operators 


{} 
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Covariant (exterior) derivative, Riemannian covariant (exterior) derivative 
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Convective covariant derivative (see, e.g., eq. (|55(l) 
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